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Abstract 



In this paper, we prove that the 5-dimensional Schwarzschild-Tangherlini 
solution of the Einstein vacuum equations is orbitally stable (in the fully non- 
linear theory) with respect to vacuum perturbations of initial data preserving 
triaxial Bianchi-IX symmetry. More generally, we prove that 5-dimensional 
vacuum spacetimes developing from suitable asymptotically flat triaxial Bianchi 
IX symmetric initial data and containing a trapped or marginally trapped ho- 
mogeneous 3-surface necessarily possess a complete null infinity T"*", whose 
past J^(X^) is bounded to the future by a regular event horizon 7i^, whose 
cross-sectional volume in turn satisfies a Penrose inequality, relating it to the 
final Bondi mass. In particular, the results of this paper give the first examples 
of vacuum black holes which are not stationary exact solutions. 

1 Introduction 

The study of higher-dimensional gravity has attracted much attention in recent 
years, motivated mainly by speculations from high energy physics. The variety of 
possible end-states for vacuum gravitational collapse in higher dimensions appears 
richer |S] than in 4 dimensions and gives rise to many interesting questions. All 
analytical work, thus far, however, has centred on the question of the existence and 
uniqueness of static or stationary ^1^] solutions, or has been based on study of 
the linearized equations ^2 El- While such results are suggestive as to what may 
occur dynamically, they do not directly address the problem of evolution and leave 
open the possibility that the non-linear theory admits phenomena of a completely 
different and unexpected nature. 

The purpose of this paper is to initiate the rigorous study of dynamical vacuum 
black holes in higher dimensions in the fully non-linear theory. Specifically, we will 
study the problem of evolution for the Einstein vacuum equations 



for asymptotically flat initial data possessing triaxial Bianchi IX symmetry. This 
model has been recently introduced by Bizon et. al. [2. They show that vac- 
uum solutions with this symmetry have two dynamic degrees of freedom, and the 
Einstein equations can be written (see P]) as a system of non-linear pde's on a 
2-dimensional Lorentzian quotient of 5-dimensional spacetime by an SU (2) action 
with 3-dimcnsional orbits. 
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The system of equations thus obtained is studied numerically in [2| , where ana- 
logues of eritieal behaviour have been discovered. Proving rigorously the kind of 
behaviour suggested by these numerics appears a formidable problem, beyond the 
scope of current techniques. Implicit in the discussion of [2], however, is the notion 
that there is an open set of initial data that leads to black hole formation. It is this 
aspect of 121 that we will formulate and rigorously prove in this paper. 

The main result is 

Theorem 1.1. Consider asymptotically flat smooth initial data {S,g,K) for the 
vacuum Einstein equations, possessing triaxial Bianchi IX symmetry. Let {Ai,g) 
denote the maximal Cauchy development, and let ir : M —i- Q denote the projection 
map to the 2- dimensional Lorentzian quotient Q. Suppose there exists an asymp- 
totically flat spacelike Cauchy surface S C Q, and a point p G S such that 7T^^{p) 
is trapped or marginally trapped, and (at least) one of the connected components 
S\ {p} containing an asymptotically flat end is such that Tr~^{q) is not outer anti- 
trapped or marginally antitrapped for any q in the component. Then Q contains a 
subset with Penrose diagram: 




Moreover, the null infinity X"*" corresponding to the above end is complete, and the 
Penrose inequality 

r < a/2M7 

holds on TC^ , where r denotes the volume-radius function, and where Mf denotes 
the final Bondi mass. 

Note that one can construct a large family of initial data such that the assumptions 
of the theorem are satisfied with S = 7r(iS), for instance. 

The region J^(X+) depicted above is what is typically called a black hole exte- 
rior, the region Q \ (X+ ) is called the black hole, and TC^ is the event horizon. 
Thus, the statement of the theorem can be paraphrased by 

Asymptotically flat triaxial Bianchi IX-symmetric spacetimes evolving from suit- 
able data, with an SU (2) -invariant trapped or marginally trapped 3-surface, possess 
a black hole with a regular event horizon (satisfying a Penrose inequality) and a 
complete null infinity. 

Theorem ll.ll can in fact be specialized to yield 

Corollary 1.1. Let {S,g,K) denote initial data evolving to the Schwarzschild- 
Tangherlini metric. Then for smooth triaxial Bianchi IX-symmetric initial data 
{S' ,g' ,K'), sufficiently close to {S,g,K) in a suitable norm, the result of the pre- 
vious theorem holds for the maximal development {M' ,g'), and moreover, the black 
hole exterior of [Ai' , g') remains close in a suitable sense to Schwarzschild-Tangherlini. 

CoroUarv II . II can be paraphrased by the statement: 

Schwarzschild- Tangherlini is orbitally stable within the class of triaxial Bianchi IX- 
symmetric spacetimes. 
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The results of this paper can be thought to complement previous results of Gib- 
bons and HartnoU suggesting linear stability^, and also to the uniqueness of 
Schwarzsehild-Tangherlini as a static black hole vacuum spaectime 0. Finally, we 
note that Theorem 1 1 . 1 1 gives in particular the first examples of vacuum black holes 
which are not static or stationary exact solutions.^ 

2 Triaxial Bianchi IX 

We will say that a globally hyperbolic spacetime {A4 , g) admits triaxial Bianchi IX 
symmetry if topologically, Ai ~ Qx SU{2), for Q a 2-dimensional manifold possibly 
with boundary, and where global coordinates u and v can be chosen on Q such that 

g - -n^du dv + (^e^^af + e^^a^ + e-2(B+C)^2^ (2) 

where B, C, O, and r are functions Q ^ M, and the are a standard basis of left 
invariant one-forms on SU(2), i.e. such that coordinates (0, 0, ^) can be chosen on 
SU{2) with 

a I = sin 6 sin ijjdcj) + cos ijjd0, 

(72 = sin 6 cos Tpdcf) — sin ipd6, (3) 
(T3 — cos 9d(j) + dip. 

If there is a boundary T to Q, it is to be a timelike curve, characterized by r 0. 

From the above, it is clear that the metric Q admits an SU{2) action by 
isometry. The boundary F corresponds to fixed points of the group action. We call 
it the centre. The angular part of the metric can be understood as a "squashed" 
3-sphere. In the ease that B ~ C, the so-called biaxial case, the system enjoys an 
additional U{1) symmetry. If i? = C = we have 50(4) symmetry and the unique 
solution to the Einstein vacuum equations is five-dimensional Schwarzsehild, which 
we will here refer to as the Schwarzsehild-Tangherlini solution. 

From the Einstein equations we derive the following equations: 

a„ {n-'d^r) = + + (CO^) , (4) 

a„ {n-^d,,r) = - + S,„C„ + (C„)2) , (5) 

-2dud,, log n ~ -r.„„ = B,, (2B,„ + C7,„) + C., (2C „ + , (6) 

r ... 

dud, logl7-K-r.„.„+3^^ = -^-l (B, . + C,J + C,„ (2C,„ -t- B „)) , (7) 

r 2r^ 2 

where p denotes the scalar curvature of the group orbit: 

= e2B+2C ^ g-2B ^ g-2C _ ig-(4B+4C) _ 1^45 _ 1^40^ (g) 

1 2 2 2 ^ 

^See also Ishibashi and Kodama 1131 . 

^Such solutions are yet to be constructed in 3 + 1-dimcnsions, as, in view of Birkhoff's theorem, 
it is impossible to reduce the problem to a 1 -|- 1-dimensional system of pde's. Solutions with a 
future complete, but not past complete, X+ have been constructed, however, by Chrusciel by 
solving a certain parabolic problem. 
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From these equations we can derive a system of nonlinear wave equations for the 
four quantities r, 57, B, and C: 

1 O^p 2r,„r,„ 



f.v.v „ J (J) 

S r r 



dud, \ogn^^ + \r,ur., - \ (B,, (2B „ + C „) + C,, (2C „ + B „)) , (10) 
2r^ 2 

-S,tii, = —T,—B,v — 7;~B,u (11) 

2 r 2 r 

oT 

3 7" 3 7" ^ 

2 r 2 r 

Note that the last two equations become identical in the biaxial case. Equations 
Q and ^ are to be thought of as constraints which are preserved by the evolution 
of ©-(112). 

The system of equations (P|- H12|l should be compared to the equations originally 
derived in [2] in r, t coordinates. 



3 The initial value problem 

Consider an asymptotically flat triaxial Bianchi IX vacuum initial data set"^ (5, g, K). 
Let {M,g) denote the maximal development of {S,g,K). By standard arguments, 
it follows that {A4 , g) is triaxial Bianchi IX symmetric in the sense of the previous 
section. Moreover, the range of the null coordinates can be chosen to be bounded, 
defining i.e. a conformal embedding of Q into a bounded subset of ]R^+^. The two 
possibilities for the global structure of the image of such an embedding are depicted 
below: 




depending on the number of asymptotically flat ends. S above denotes Tr{S). In 
what follows, the notations J+, closure, etc., will refer to the topology and causal 
structure of M^+^. By the definition of asymptotic flatness, it follows that r tends 
monotonically to infinity along S, sufficiently close to the points labeled iq. More- 
over, Qn J^{S) is foliated by constant-w curves emanating from S, and constant-u 
curves emanating from 5* U F. 

^We leave to the reader the correct formulation of this notion. 
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4 Local existence and extension 

We wish to understand those TIPs in Q which do not "emanate" from the centre 
r. For this, the foUowing local existence theorem in null coordinates shall suffice 
for our purposes. 

Proposition 4.1. Let ^l, r, B, and C be functions defined on X = [0,d] x {0} U 
{0} X [0,d]. Let k>0, and assume r > is C^+'^{u) on [0, d] x {0} and C'=+2(d) 
on {0} X [0, d], assume that Q,, B and C are C'^^^ {u) on [0,d] x {0} and C'^^^{v) 
on {0} X [0, d]. Suppose that equations hold initially on [0, d] x {0} and 

{0} X [0,(i], respectively. Let denote the C^^{u) norm of f on [0, d] x {0}, 

\f\n,v the C"(w) norm of f on {0} x [0,d], etc. Define 

N = sup{|r!|i,„, |r|2,„, |r|2,„, Ir^^lo, |C|i,„, 

Then there exists a 5, depending only on N, and a C'^"''^ function (unique among 
functions) r and C*'"'"^ functions (unique among functions) fl, B, and C, 
satisfying equations ^~ ^1S^ in [0,(5*] x [0,(5*], where S* = min{d, (5}, such that the 
restriction of these functions to [0, d] x {0} U {0} x [0, d] is as prescribed. 

Proof. The proof is by standard methods and is omitted. □ 

From the above Proposition and the maximality of the Cauchy development, 
the following extension principle follows. Given a subset Y C Q \ P, define 

N{Y) = snp{\nU, \n-'\o, |r|2, |r"i|o, |C|i}, 

where, for / defined on Q+, denotes the restriction of the C'^ norm to Y. 

Proposition 4.2. Let p e'Q\T, and q e Qr\I^{p) such that J^{p)nJ+{q)\{p} C 
Q, and N{J~{p) n J+{q) \{p}) <oo. Then p e Q. 



5 The Hawking mass 



A remarkable feature of the system of equations l(^- H12() is the existence of energy 
estimates for B and C. For this, we first define the so-called Hawking mass 
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We compute the identities: 
4 



du-m 



3 02 



4 r 



3 02 

Note that p is bounded above: 



2 

1 - -p 



2 ' 



3 



(13) 

(14) 
(15) 
(16) 



(A straightforward way to show this is to set x = e'^^ , y — e^'-^ and to study the 
function p{x,y). First one shows that p{x,y) < | in the region 

r 1 11 

R 



X < -^,y < a {x>10,y> 10}. 



10' 



10 



(17) 
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Next one determines the critical points of p{x,y). It turns out that that there is 
only one extremum at x = 1, y = 1, which is shown to be a maximum. This proves 
p{x, y) ^ f with equality only for the round sphere, B = C = 0.) 

By ()16|l . we now see that all terms in square brackets are manifestly non-negative. 
Thus, if, say dur < and dyr > 0, we have 

dum < 0, a„TO > 0. (18) 

The monotonicity (|18|l can be compared with the monotonicity in the r-direction 
for the (r, t) coordinates given by formula (8) of 

6 The regions 7Z, T, and A 

Let us define the regular region 

n = {pe Q such that dyr > 0, a„r < 0} (19) 

the trapped region 

T = {peQ such that dyr < 0, (9„r < 0} (20) 

and the marginally trapped region 

A^{pe Q such that d^r = 0, d^r < 0}. (21) 

The reader is warned that the term regular is meant with reference to the asymp- 
totically flat end in the direction of which the vector dy points. By the results of 
the previous section, the inequalities (|18|l hold in TZU A. In the next section, we 
will show how this leads to a stronger extension theorem than Proposition 

7 Extension in the non-trapped region 

The monotonicity (|18|l indicates that our system l|H)- (|12|) shares a formal similarity 
with spherically symmetric 3 + 1-dimensional Einstein-matter systems, for suitable 
matter fields satisfying the dominant energy condition (See [SJE]). In particular, 
one might conjecture that an extension principle analogous to the one formulated 
in [2] holds in the non-trapped region. This is what we show in this section. 
We have 

Proposition 7.1. Let p eTZ\T, and q eTZUAnl^{p) such that J"(p)n J+(g)\ 
{p} C 7^ U A. Then peTZUA. 

Proof. The proof adapts techniques introduced in [5] . Let us introduce the following 
notation: 



1/ = dyr, 


(22) 


A = dyr, 


(23) 




(24) 


4 ly 


2m 


(25) 




Cb = riduB, 


(26) 


Cc = ridyC, 


(27) 


9b = ridyB, 


(28) 


Oc = ridyC. 


(29) 
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Note that k{1 — /i) = A. 

Let X denote \ I~^{q)) H Q. Setting p = {ui,vi), q = [u^.v^), we have 

X = {ite} X [iiciii] U [uejiti] X {I'e}. Since X is compact, the quantities 

r, K, A, V, m, B, C, CbXc, &b, Oc, du^, dM, 9„A, (30) 

are uniformly bounded above and below on X: 

0<rQ<r< R, 
< A < A, 
> iyo> v > N, 

\B\<Pb, 

\C\ < Pc, 

\eB\ <Tb, 

\0c\<Tc, 

\Cb\ < Zb, (31) 
\Cc\ < Zc, 
\m\ < M, 

< Ko < K < K, 

\dun\ < H, 
\d^n\ < H, 

\duiy\ < H, 
\dyX\ < H. 

By Proposition 14.21 in view also of the fact that il^ — ~Anv, to prove Proposi- 
tion it suffices to show that the quantities (|30() are uniformly bounded every- 
where in [mejUi] X [ue,t;i] \ {(uijfi)}, with bounds similar to (|31|l . 

We first derive a bound for r. Integrating v along u, and A along we obtain 
from ipT)l . in view of the signs of v, A in 7^ U A, that 

< ro < r < i? (32) 

in X [i'e,wi] \ {(mi,wi)}. a similar argument can be given for the mass: 

Integrating H14|) along u yields 

m{u*,v*)-m{u,,v*)<Q (33) 

and integrating (|15|l along v yields 

m{u,v*)-m{u*,v^)>Q. (34) 

We conclude the bound 

-M <m<M (35) 

on [u^,ui] X [ve,vi] \ {(ui,'!;i)}. 

A bound on k can be derived as follows: Note that k > by definition, in view 
of the v <0. On the other hand, we compute from Q 

«^,u = -l^^'' {{B,uf + B^uC.u + {C.uf) < 0. (36) 
Thus, integrating in u from X, in view of H31|l . we obtain 

< K < A' (37) 
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in [we, ui] X [ve,vi] \ {{ui,vi)}. 

Next wc bound the quantity v using the evolution equation (Q, written: 



1 n^p 

3~ 



(^v^ — ^ ,UV ' 

Integrating this equation in v we get 



2iyX _ Ukp 2X\ 
r \ 3r r J 



/^(u ,v ) = i^[u ,Ve)exp \J [— — \ (u ,v)dv 



Since p < |, A > 0, we obtain the upper bound 



-I' < \N\ ■ exp 



2eK 



= N'. 



(38) 



(39) 



(40) 



From the above and (|37|) it follows that the quantity fi^ = —Akv is also bounded 
from above. 

To estimate B and C, we revisit the equations p4|) and \ib\ . in view of H35|l . 
to infer a-priori integral estimates for derivatives of these quantities. Equation (|15|) 
gives 

'^i^-^ {{B.^f + + (C„)') +\rU- {u\ v) dv < 2M (41) 



and therefore, since 



BI + ByCy + Cl> -BI + -Cl > 0, 



we have 



3 K 



[B^^f {u*,v) dv 



3k 



;) {u*,v) dv < AM. 



(42) 



Obviously, the same inequality holds with B replaced by C. In the same way, 
integrating equation H14|l along u using the mass-bound H35() leads to the estimate 



(43) 



Again, the same inequality holds with B replaced by C. 

We may now integrate the equation B^y — r~^9B in v to obtain 



\B{u\v*)\ < \Biu\v,)\ + 



<Pe 



— (u* , v) dv 



^dvJ I :Ldv<PB+VT2Mj^ = Pi 



(44) 



where we used the Schwarz inequality in the step from the first to the second line and 
(|42|1 for the last step. In a completely analogous fashion-integrating Ct, = r^^Oc 
in w-we obtain the same bound for C . Having bounded B and C, it follows from 
© that p is also bounded in [ug, ui] x [v^, vi] \ {(ui, ui)}. This enables us to bound 
A. Rewriting the evolution equation Q for r^uv in terms of quantities we already 
control we obtain 



9„A = r. 



An 2 , 
— p- -K I 1 
6r r 



2m 



(45) 
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which we can integrate along u. Because we already control all the quantities 
appearing in the integrand we immediately obtain a bound for A in [u^, ui] x [v^, 

\<L. (46) 

The determination of a suitable constant L is left to the reader. 

We turn to bound and n from below, away from zero. In view of the bound 
on IpI, we may derive immediately from 1)39(1 a bound 



For K, we integrate rewritten as 

dun = n {Kv-\Bl + B^Cu + Cl) 



(47) 



to obtain 



u 2 



> 



where, for the last inequality, we use and the bounds proved above, in partic- 
ular, the M-analogue of l|lT)) . 

Finally, we note at this stage that from it follows immediately r^uv is 
bounded in [u^, Ui] x [v^, vi] \ {(ui, ■^i)}. 

We turn now to bound the derivatives of B and C. First let us consider d^B, 
dyC: Differentiating Ob = r^dyB in u and using the evolution equation Hll|) we get 



2 r 3vr 



(48) 



which can be integrated in u to give 



\eB{u\v*)\ < \9Biu,,V*)\ + ^ 



I — ^ (u, t;*) du 



3Vf 



2B+2C , „-4S-4C 



2e 



(49) 



The term in the last line is bounded because we control all quantities in the in- 
tegrand. We estimate it say by the constant F. For the second term we use the 
Schwarz inequality and the a-priori bound (|43|l : 



\0b{u\v*)\<Tb + F+- 



(m, v"^) du 



<Tb+F 



(_^)(1_^)(^) (^u,v-)du 



(u, V*) du 



<Tb + F+ ^VuM- K ■ \lr^^ + Mr^'^ = V. 

(50) 



Hence we bounded 6b and therefore dyB. The bound for 9„C is obtained completely 
analogously. 
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Next we turn to duB, 9„C. Differentiating (b = r^duB with respect to v using 
the evolution equation we obtain 

It o y 7* 

Integration in w now yields a bound for since all the quantities on the right 
have already been shown to be bounded. (Alternatively we could use the Schwarz 
inequality and the a-priori bound (|42|) .') The bound for Qc and therefore C„ is 
obtained in a completely analogous manner. Having bounded B, C and their first 
derivatives, equation yields that B,uv is also bounded. 

Bounds for r2„ and ^l^ follow by integrating (|10|l in v and u respectively. Finally, 
bounds for r.„„ and r_^y follow from Q respectively |SJ) and the previous bounds. 

As remarked at the beginning, the proof now follows by applying Proposition l4.2l 

□ 



8 Null infinity 

Let S be as in the statement of Theorem ll.il Without loss of generality, let the 
asymptotically flat end in question be such that 9„ points "outwards" . We define a 
set X+ C (Q\ Q) n J+(5), as follows: Let 

U = lu : sup r{u,v) = oo > . (52) 

[ V ■ («,t))GQ+ J 

For each u U, there is a unique v*(u) such that 

iu,v''iu)) e(Q\Q)nJ+{S). (53) 

Let the end in question have limit point on S given by io — {u.,V). Then the 
null- infinity corresponding to iq is defined as the set 

U {u,v*{u)). (54) 

uGU:v*(u)=V 

Standard arguments show that X+ is non-empty for the data considered here. It is 
straightforward to show, adapting , that X+ is then a connected ingoing null-ray 
with past-limit point io. Denote the future limit point of by . A priori, it 
could be that i+ e X+ . 

Adapting one shows from H18|) that the Bondi mass 

M{u) = limsupM(u,?;) 

is a finite (not necessarily continuous) function on non-increasing in u. We 
define Mf = inf M[u) to be the final Bondi mass. 

9 Proof of Theorem 11.11 

This proof is an adaptation of methods introduced in 

As above, let S be as in the statement of Tlieorem ll.il let dy be the outward 
direction, and consider the set 

V = J+ {S) c^ J- {!+) f^ Q. 
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This set is non-empty. On the other hand, by the Raychaudhuri equations Q-©, 
and the assumption that 9„r < along S, it follows that 9„r < along future- 
directed constant-w curves in Q emanating from SO {v > v{p)}, and 

VcTZ. 



Since by assumption p & T U A, it follows that p ^ V, and thus V has a non-empty 
future boundary in Q. Denote this boundary 7^+. Note also that m > r^{p) > in 
V, and thus in particular, Mf > 0. 

Proposition 17.11 shows immediately that Ti.'^ cannot terminate before reaching 

i.e., the Penrose diagram is as: 




We will first show that the latter is the case, i.e. z+ ^ X+, in fact, that the Penrose 
inequality 

<2Mf, (55) 

holds on the event horizon 7i+. 

To show (|55f) on one assumes the contrary, i.e. the existence of a point 
iU,V) with r'^{U,V) = > 2Mf on the horizon, and as in one infers (using 
monotonicity properties of r and m, together with ProDOsition l7.1|l the existence of 
a neighbourhood of the horizon which is part of the regular region: 

Af -.^ [ua,u"]x [V,V) an (56) 

with uo < U < u" . In particular this neighbourhood can be chosen such that there 
exists an R' < R with the property that in U, u" x [1/, V) <ZTi 

2m '^M 

r> R' and 1 ^ > 1 - (57) 

r-2 - (i?')2 

holds. The last step is to show that for any u* £ [uq^u"], limy*^rx>'r{u* ,v*) — oo, 
i.e. Ti"*" cannot be the event horizon, as defined, a contradiction. 

To show this last step, having shown (|57|l . we proceed as follows: Integrating 
(|14|l along u from uq to a point u* < u" we obtain the estimate 

v>V •'^0 



sup / ((B„)2 + B^uC,u + (C„)2) {u,v) dH < M (58) 



which can be written as 



sup/ t(1"M)((^,«)'+B„C.„ + (C„)2)(ii,T7) cRl< A/. (59) 

Taking H57|) into account we can derive the estimate 

sup/ -— -{u,v)cM>-—^ — — (60 
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valid for any u* G [uo,u"). Integrating (|47|l . wc obtain 



—2M 

k{u\v*)>k{u^,V*)-C1,Y'[ ) (61) 



and therefore 

Ku\ > ^1 - ^ j exp ^ j A(.o, .*)• (62) 

Integrating (|62|l in v, we see that hmt,*_»v r(u*, w*) ^ oo, since hnit,*^y ?'(wo, f*) ^ 
oo on the right by the definition of X+. We conclude (it*, V) £ I'^ . Therefore, Ti"*" 
is not the event horizon and we have arrived at the desired contradiction. 

The only thing left in the proof of Theorem II. II is to show the the completeness 
of X+. (Completeness here refers to an adaptation in [S] of the concept defined 
in [3].) We have to show that the suitable normalized affine length, as measured 
from a fixed outgoing null curve u = tig, of the ingoing null-curves v = const in 
J^(X+) tends to infinity as v . More precisely, we define the vector field 

iVyu.v) au 

on J^(X+) n Q+. Note that this vector field is parallel along all ingoing null-rays 
and along the curve u = uq. We will show 

rU _^ 
lim / {X{u,v)-u) du = oo. (64) 

From equation Q we can derive 

r2^(u, t')r2~^(Mo, f) — i'{u,v)i'~^{uo,v) ■ (65) 



■ exp ^ 1^ ((S,„)2 + B „C„ + (C„)2) {u, v)du 



Let M be the Bondi-mass at uq. We choose an R such that R^ > 2M > 2Mf and 
consider the curve {r = R} C] (X+ ) . For sufficiently large vq < V, all ingoing 
null-curves with v > vq intersect {r ~ R} fl J~(I+) at a unique point {u*{v),v), 
depending on v. 

Analogously to H6U|) we derive the bound 

" 2r —9M 
- {{B.^ f + B,,Cu + (C „)2) {u, v)dM > , (66) 

Uq 

which we use to estimate 

U pu* {v) 

{X{u,v)-u) du > {X{u,v)-u) du 



Uo 



= v^'^{ua,v) I exp 



r(Mo,w)-i? / -2M 



{-v){ua,v) '^\R^-2M 



(67) 



Since r(uo, w) — > oo as w — > cxd we only need to show that {—v)[uq, v) is uniformly 
bounded in v. The quantity 
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satisfies 

dvj^ = + + (69) 

1 — f-i 1 — fi OA 



which integrates to 



exp 



1 — /i 



y i /I 



We can choose vq (so large) such that 

2M 

{r{uo,VQ)) 

Set i?' = r(Mo,wo)- Analogously to (|60() and H66|l we derive the bound 



which enables us to obtain from 170() the estimate 



(70) 



, , 2i\/ \ ^ / 2M , 
-K"o, «) < 1 - TT^TTI exp (73) 



(i?')V V(^')^-2Af 

for w > 1^0, which in turn shows uniform boundedness of (—i'){uo,v) in v. 

We have shown above the future completeness of T~^. The past completeness 
follows by standard results assuming sufficiently tame initial asymptotic behaviour. 
In particular, it is immediate in the case that B, C, VB, VC vanish initially outside 
a compact set. 



10 Proof of Corollary [TT] 

Let 5* denote the projection of an arbitrary spherically symmetric Cauchy surface in 
Schwarzschild, and let S denote the projection of a second asymptotically ffat spher- 
ically symmetric Cauchy surface, with the property that 5* contains a p satisfying 
the conditions of Theorem ll.il (Such Cauchy surfaces clearly exist.) By Cauchy sta- 
bility, sufficiently small triaxial Bianchi IX perturbations of Schwarzschild data on 
TT~^{S) yield solutions {A4',g') possessing a triaxial Bianchi IX symmetric Cauchy 
surface S' with geometry arbitrarily close to that of 5, in particular, also satisfying 
the assumptions of Theorem ll.il We apply thus this theorem. 

Finally, we note that the Hawking mass on 5" is arbitrarily close to the constant 
value M it takes on Schwarzschild, i.e. we have M — e < m < M + e on S. By 
the monotonicity (|18|l . it follows that this bound is preserved in J+(S") n J^(X+). 
It is this statement-together with the stability of the Penrose diagram and the 
completeness of X+ -that we term "orbital stability". 



11 Final comments 

Besides orbital stability, one is interested in what could be called asymptotic stability 
of the Schwarzschild family, i.e. the statement that perturbations of a Schwarzschild 
initial data set asymptotically approach another Schwarzschild solution. An even 
more ambitious problem would be to understand the rates of approach, as in |7|. 
These problems remain open. 
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Another interesting and partly related question is to understand the structure 
of the outermost apparent horizon.^ In analogy to (S], we may define this as the set 



A' = {{u,v) e A : iu*,v) e 7^for all u* <uand3u' : {u',v) G J~ {I+)r\QnJ+ (S)} 

(74) 

As in [S], A' is now easily shown to be an achronal curve intersecting all ingoing 
null curves for v > vq for sufficiently large vq. In addition, one shows easily that on 
A' , the Penrose inequality (|55() holds. There are many other issues, however, which 
are not settled: Is it a connected set in a neighborhood of i+? Is it "gcnerically" 
a strictly spacelike curve in a neighborhood of i"*"? Does it terminate at i+ in the 
topology of the Penrose diagram? For more on these questions, the reader can 
consult the literature on so-called dynamical horizons, in particular pQ. 
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